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Abstract 

An extension of potential theory in R" is obtained by continuing the Euclidean distance function holo- 
morphically to C". The resulting Newtonian potential is generated by an extended source distribution 
S{z) in C" whose restriction to R" is the point source S{x). This provides a possible model for extended 
particles in physics. In C"^^, interpreted as complex spacetime, 5 acts as a propagator generating so- 
lutions of the wave equation from their initial values. This gives a new connection between elliptic and 
hyperbolic equations that does not assume analyticity of the Cauchy data. Generalized to Clifford anal- 
ysis, it induces a similar connection between solutions of elliptic and hyperbolic Dirac equations. There 
is a natural application to the time-dependent, inhomogeneous Dirac and Maxwell equations, and the 
'electromagnetic wavelets' introduced previously are an example. 



1 Motivation and Preliminaries 

Most fundamental theories of physics are based on the concept of potentials and fields generated by point 
sources, which presupposes that objects or "particles" can, in principle, be localized within arbitrarily small 
regions of space and/or time. This is a vast extrapolation from empirical evidence, and it should perhaps 
not come as a surprise if such theories experience some fundniental difficulties. In Newtonian mechanics, 
the problem of N point "bodies" interacting through gravitation has, in general, no solution due to the 
possibility of collisions. This becomes a serious difficulty for iV > 3, where the set of initial conditions 



leading to collisions is nontrivial |AM78, Chapter 10]. In classical electrodynamics, difficulties arise where 
the field produced by a point charge unavoidably acts back on the same charge, leading to infinite self- 
energies and run-away particle trajectories [ J99|. A way out of this dilemma was proposed by Wheeler and 
Feynman [WF4£, WF49|, but their action-at-a distance theory, apart from being highly counter- intuitive 
[F64, p. 28-8], has resisted quantization and is not generally regarded as being a fundamental description 
of Nature. In quantum electrodynamics and other quantum field theories, point particles cause divergences 
which necessitate infinite "renormalization" procedures, a subject of some contraversy [ C99{ . String theory 
jP98| does, in fact, not need infinite renormalization because its basic objects (strings) are extended in space 
rather than mathematical points. This is one of the reasons it is regarded with great hope as a possibility 
for unifying physical theories. However, a full development of (super)string physics is rather difficult and 
not expected to near completion for many years. 

The ideas developed here began many years ago, motivated in part by the hope that an extension of physics 
to complex spacetime, justified at the foundational level, might give a way to circumvent the problems 
associated with point sources by applying residue methods. After some years of study and research this 
led to papers IotI ^7|, ^80[ and books p?90| , ^9| whose main thrust has been to develop a direct 

physical interpretation of the complex spacetime as an extended phase space, with the imaginary spacetime 
parameters carrying directional information while the real spacetime parameters describe (approximate!) 
localization. A major stumbling block in this program has been the construction of extended sources, since 
it seemed that they tend to spoil the holomorphy of the theory globally rather than just locally. Here we 
propose a natural solution to this problem. It turns out that although holomorphy enters the theory at the 
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level of the fundamental potential or Green's function and this indeed gives a canonical extension of general 
fields to complex spacetime through convolution, these fields need not be holomorphic anywhere due to the 
specific structure of the extended sources. 

We begin in Section ^ by extending the Euclidean distance r{x) to a function j{x + iy) holomorphic in a 
domain On C C". This leads to a natural coordinate system in R" which will play an important role in the 
sequel: the oblate spheroidal coordinates adapted to y ^ 0. 

In Sections ^ and ^ we use j{z) to extend the Newtonian potential (j){x) (fundamental solution of the 
Laplacian in K") to On- The resulting holomorphic potential (f>{z) has a source distribution S{z) in C" 
which, although nowhere holomorphic, is nevertheless a canonical extension to C" of the point source S{x) 
in K". For even n > 4, S{x + iy) with fixed y ^ is supported in x on the sphere B{y) of codimension 2 
and radius |y|, centered at the origin and lying in the hyperplane orthogonal to y. For all other values of n, 
it is supported in a; on a membrane whose boundary is B{y). The membrane is determined by a choice of 
branch cut in 7, and in the simplest case it is a disk Ea{y). 

In Section^ we compute the sources S{x + iy) explicitly as distributions in and supported, respectively, 
in Ea{y) and B{y). This is especially interesting in the physical case of R"^, where 5{x + iy) consists of a 
uniform line charge on the rim of the disk Ea{y) accompanied by a simple layer (with zero net charge) and a 
double layer distributed on the interior. Hence the original unit "charge" carried by 6{x) is now spread over 
a circle of radius and a variable polarization density is induced on the disk swept out by "blowing up" 
the point source from the origin to its rim. That so much natural structure results from a simple analytic 
continuation is remarkable. 

In Section ||we show how to include time in this framework. Since we already have a holomorphic extension 
4>{x, s + it) of the fundamental solution (t){x, s) of the Laplacian in R""*"^, it is natural to interpret i as a time 
parameter and look for a connection with the wave operator A3; — d^, which is the "analytic continuation" 
to R"'^ of the Laplacian A^, + in R"+^, now interpreted as Euclidean spacetime. We show that d{x, s + it) 
acts as a propagator, generating solutions of the initial-value problem for the wave equation in R"'^ from 
Cauchy data at t = 0. This gives a connection between the Laplace equation and the wave equation without 
assuming analyticity of the initial data, because the propagator S{x, s + it) is not holomorphic in s + it. 

In Section |^ we extend our method to Clifford analysis, where the Laplacian and the wave operator are 
replaced by their "square roots," the elliptic and hyperbolic Dirac operators. The time-dependent Maxwell 
equations are a natural application of the ensuing formalism. 

The idea of extending harmonic functions to obtain solutions of the wave equation has been studied by 



Garabedian |G64| and, in a more general context, by Ryan |R90, R90a, R96a|. However, the methods 
employed in these references depend critically on the assumption that the boundary data for the harmonic 
functions is holomorphic. As emphasized above, our method is free of this restriction. 



2 Complex Distance and Spheroidal Coordinates 

For n>2, let 

z = x + iyeC'', \x\=r, \y\ = a. (1) 

Define the complex length of z as 

7(2;) = Vz^ = \/ x'^ — y'^ + 2ix ■ y = \/ r'^ ~ a'^ + 2ix ■ y = p + iq. (2) 

The complex distance between any two points zi, zi G C" is then defined as 7(21 — Z2). To complete this 
definition, a branch of the complex root must be chosen. The branch points form the null cone 

M ={z \ 7(2) = 0} = {cc + iy I = y2 and a; • y = 0} , (3) 

a manifold of real dimension 2rt — 2 in C" « R^". In the context of the wave equation, M will be seen to be 
related to the light cone. We will be interested in the restriction of A/" to a fixed nonzero y e R", i.e., 

B = B{y) = {xeW\r^a, x-y^O}. (4) 
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B is the sphere of dimension n — 2 in R" obtained by intersecting the sphere r = a with the hyperplane 
orthogonal to y. We call B{y) the branching sphere with axis vector y. 

We now concentrate on j{x + iy) as a function of x, regarding y as a fixed nonzero vector. In order to 
make this function single- valued, it is necessary to introduce a branch cut in M". To see how this can be 
done, fix any x G B. In the plane determined by the two orthogonal vectors x and y, draw a circle of radius 
e < a centered at x. Each time we go around this circle, 7 changes sign. In order to obtain a single-valued 
function for 7, it is therefore necessary to cut every such circle. Furthermore, this must be done subject to 
the requirement that 7 must be an extension to C" of the usual length function in M", i.e., 

y ^ ^ 7(0; + iy) r{x) > 0. (5) 

The simplest cut (but ccrtinainly not the only one) is obtained by requiring that 

p EE Re 7 > 0, (6) 

which means that each of the above circles is cut at the point 

ajc = (a — e)x, x = x/r. 

The totality of such points, as x varies over B and < £ < a, together with B itself (e = 0), forms the set 

E„ = E,{y) = {xeW'\r<a,x-y^O}, (7) 

which is the disk of dimension n — 1 in M" obtained by intersecting the ball r < a with the hyperplane 
through the origin orthogonal to y. The boundary or rim of E„ is B: 

dE„ = B. 

(For n = 3, i?o is indeed the disk of radius a orthogonal to y and B is its rim.) In the context of holomorphic 
potential theory. Eg or B, depending on whether n is odd or even, will represent the support of a source 
distribution. As y ^ 0, both sets contract to the origin, the support of the usual point source. 

We now derive the basic properties of the real and imaginary parts p, g of 7, which will be important in our 
study of the holomorphic potentials. From (^) it follows that 

p2-g2^r2-a2, pq^x-y = aC, (8) 

where C is the projection of x onto y ^ 0. Define the cylindrical coordinate p by 

p =r -C =a +p -q ^ = ^ . 9 

This shows that the imaginary part of 7 is bounded by a = |y|: 

-a<q<a. (10) 

We will use {p, as part of a cylindrical coordinate system in R" with y as its "z-axis." Equations (^) and 
(^) show that the surface Ep = {x \ p = constant} with a given value of p > is an oblate spheroid given by 

^P- ^^ + \ = 1- (11) 

Similarly, the level surface {x \ q"^ = constant} with a given value of < < is a hyperboloid of one 
sheet given by 

Hq-- 2 2^1- (1^) 

a'^ — q'^ q'^ 
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Note that we have so far avoided the surfaces with p — and q = 0, ±a. These give degenerate forms of Ep 
and Hq. In fact, as p ^ +0, the oblate spheroid Ep contracts to the disk (Q): 

p^O ^ Ep^{x\(: = 0, p<a} = Eo. (13) 

More precisely, since the interior of Eq is covered twice, we will distinguish between its front and back sides: 

E„=E+UE-UB, (14) 

where 

£;± = {cc I C = ±0, p<a} = {x\p = 0, ±q> 0} (15) 

are the interiors of the front and back sides of Eg. Although i?+ and Eg coincide as sets, the distinction 
between them will be very important for the following reason. It can be easily seen that 



xeE^ => 9 = ±Va2^^= ±v42~~^. (16) 

This shows that while p is continuous across E„, q has a spherical jump discontinuity there. In the context 
of holomorphic potential theory, i?+ and E^ will be seen to form the two sides of a double layer. 

The degenerate values of q are q = 0, ±a. As g — > ±a, the semi-hyperboloids with ±g > collapse to 
the positive and negative ^-axis, respectively. As q ^ 0, Hq collapses to the set Hq = {a; | ^ = 0, p > a} 
which, like Eg, is covered twice, once by the limit of each semi-hyperboloid H^ as g ^ ±0. But in this case 
we do not distinguish between the two copies because 7 is continuous across Hq. 

As mentioned above, the parameters C and p are part of a cylindrical coordinate system in M". Note that 
the intersection of Ep and Hq is 

= E^ n Hi 



'p n iJ, = {a; e M" I aC = pg, ap = Va^+pVa^-?^} 



Here and henceforth, we denote the pair of real variables {p, q) by the single complex variable j — p + iq for 
convenience. If g 7^ ±a, then S-y is a sphere of dimension n — 2 and radius p. The variable point on S-y may 
be represented as 

pcT + Cy e Sj, (17) 

where cr runs over the unit sphere of dimension n — 2 in the hyperplanc orthogonal to y. In particular, 7 = 
gives C = and p = a, so that 

So^{a(T\(TeS''-^} = B{y). 

If n = 2, S.y consists of just two points. If n = 3, Sj is a circle. Thus, for fixed y 7^ 0, a complete set of 
cylindrical coordinates in M" is given by (p, C, cr). Since C and p are functions of (p, q), we may equivalently 
use the coordinates 

a; (p,q,cr) = (7,0-). (18) 

This is the oblate spheroidal coordinate system in M", long known as a useful tool in the theory of special 
functions [T96| . It is remarkable that this system is so intimately related to the holomorphic extension of 
the distance function. 

Now let V be the gradient with respect to x with y held constant, and let A = be the Laplacian in M". 
Then 

7^=2^ ^ 7V7 = 2;, hence (V7)^ = V7 • V7 = 1. (19) 
Taking the divergence of (|l^ ) gives 

n — 1 

V7 • V7 + 7A7 = n, ^ A7= . (20) 
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Now 



therefore 



„ z 72; „ px + qy „ py — qx 



V7-V7 = l ^ {Vpf -{Vqf = 1, Vp-Vq^O 



V7 • V7 = ^ => {VpY + [VqY = 



which gives 



(Vp) 



77 



2 _ a2 + p2 



77 



77 



9 ? 

2 ^ a - g 

77 



Finally, the real and imaginary parts of (]20|) give 



n - 1 
77 



■p, Aq 



n — 1 
77 



(21) 

(22) 
(23) 

(24) 

(25) 



We will need to compute volume integrals in the oblate spheroidal coordinates. Note that the area of the 
unit sphere S'"'^ C K" is pfe^ 



rW2)' 



n > 2. 



(26) 



Let dcr denote the surface measure on the unit sphere S"" ^ in the hyperplane orthogonal to y 7^ 0, normalized 
so that 

(27) 



da- = 1. 

Then for fixed y, the volume measure dx in M" is given in cylindrical coordinates by 

dx = LOn-i p^^^dpdQdcr. 

By (1) and (|), 



(28) 



ri-3 n-3 



a V" {pdp - qdq){pdq + qdp) = a p" '\p + q )dpdq, 



where dpdq denotes the antisymmetric exterior product of differential forms (see |GS64], for example). 
Therefore the volume measure in the oblate spheroidal coordinates is given by 



dx = p"--^ (p^ +q^)dp dq dcT = dp dq da. 

a a 



(29) 



3 Holomorphic Potentials and Their Sources 



For simplicity, we now assume that n > 3. The case n = 2 is similar but requires some special attention and 
will be described elsewhere. Consider the fundamental solution (p{x) of Laplace's equation defined by 



A(j){x) = 6{x), lim (j){x) = 0. 



(30) 



It is given uniquely by |CH62 



2-n 



1 r 

ujn 2 — n 



X e 



n > 3. 



(31) 
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For n = 3, 0(a;) — — l/Anr is the Newton-Coulomb potential with unit mass or charge. Define the holomor- 
phic potential in C" by 

cj,{z)^ — ^, ^eC", n>3. (32) 
ojn 2 — n 

For odd n, (t>{z) inherits the branch cut of ^(z). For even n, the only singularities occur on B, where 7 = 0. 
Thus 7(2;) and (f>{z) are analytic continuations of r{x) and <j){x) to the domains 

On = {zeC' \ p>0} = {x + iy\x(^ Eo{y)} , odd n > 3 

= e C" I 7 7^ 0} = {cc + I a; ^ S(y)} , even n > 4. 



Proposition 1 For fixed y, (f>(x + iy) is harmonic with respect to x when x + iy (z On- 
Proof: By (^2|), we have 

c^„V0(^) =7i-"V7. 

Thus by (|^), 

w„A(/)(^) = (l-n)7-"+7i-"A7 = 0. I 

Our objective is to compute the source distribution of 4>i^): which we define formally by analogy with ( |30| ) 
as 

J(2) = A(f>{z), z e C". (33) 



This will be shown to be a generalized function |GS64| of x for any fixed y, meaning that given any sufficienlty 
smooth "test" function f{x), the integral 

{lf)= f 5{x + iy)f{x)dx (34) 

defines a bounded linear functional of /. (Generalized functions are more commonly known as distributions 
[Z65|. We use the former term here in order to avoid confusion with the term "source distribution.") 

By Proposition |l|, S{x + iy) is supported on a; e S for even n > 4, and on a; e E'o otherwise. In any case, 
it has compact support in the variable x, and as y — s- 0, this support contracts to the origin. We will show 
that 

y ^ ^ (J,/) ^ /(O), hence ~5{x + iy) 5{x). 
To compute the generalized function 5{z), we first define the regularized potential: 

Mz)=e{p-e)cf>{z), e>0, where 0{O = T ^ ^ (35) 

10 it 4 < 

is the (Heaviside) unit step function. The regularization ( p5| ) eliminates the singularities on Eg and B, 
replacing them by a discontinuity in (j>^ across the spheroid E^. Thus, while the source of the singular 
potential (t>{z) is concentrated on Eq or B, that of the regularized potential (j)e{z) is concentrated on E^. 
The advantage gained is that while (t){x + iy) is infinite for x ^ B C Eq, (l)^{x + iy) remains bounded in a 
neighborhood of E^. 

Now 4>£{z) vanishes in the interior of E^ but is identical to (l){z) in the exterior. Therefore its source 
distribution, defined as a generalized function by 

~6,iz) = AMz), (36) 
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represents an equivalent source distribution on x G whose potential field simulates that of S{z) in the 
exterior of but vanishes in the interior. We call the regularized source distribution, and will define 

the singular source distribution 6{z) as the limit oi 6e{z) in the sense of generalized functions, for any fixed 
y e M". That is, 

ilf) = lim (4,/) (37) 

for every test function f{x) in M". Using 9'{£_) — S{£_), (^) gives 

V(/)e = S{p - e) Vp + e{p - e)V(j), 

hence 

A0e = S'{p~e)^{Vp)^ +2d{p~e)V(l)-\7p + S{p-e)(t)Ap + e{p-e)A^. (38) 

By Proposition |l|, vanishes for p > e/2, so the last term in ( ^ ) vanishes identically. The remaining 
terms show that S^{x + iy) is indeed supported on Eg^ as expected. Inserting 

V0 = (/''(7)V7 = ^ 



into ( p3q ) and using (p2[)-(25), we have 

d"^ 4- rj-^ fi — t 

(5e(2) = [S'{p-e)(f> + 2S{p-e)(t)'] — + 5{p-e)(f> p. (39) 

77 77 

This expression represents a generalized function of a;. To make sense of it we must apply it to a test function 
f{x) in K", assumed to be sufficiently smooth. (That is, / is assumed to possess all derivatives which the 
ensuing computation requires it to possess. As will be seen, the required degree of smoothness increases with 
n.) For a given value of y G K", acts on / as in (^), 

(4,/) - / l{x + ty)f{x)dx. (40) 



Using the oblate spheroidal coordinates (|18|), let us write 

fix) = /(p<T + C y) - q, t) = /«(7, <T), 

where the two expressions on the right are obtained by substituting (^) and (|^) for ^ and p in terms of (p, q). 
Let 

Ph)= [ ni.cT)da^ [ fipcT + Cy)da^f{p,0. (41) 

The notations /''(7, cr) and /•(7) = /"(P; ?) a-re used for convenience and are not meant to imply analyticity 
in 7. Because of the normalization (^3); ./"(7) ^-nd f{p,C) are the means of fj^, er) and /(per + Cy) over 
the sphere S-y = Ep H Hq — {x : \x — (cr\ = p}. Using the expression (|2^ ) for da;, (|40| ) becomes 

/•OO /"Q 

ilj)^"^ dp dqp--^^^Ul)I\l), 



where we have used the fact that 6^{z) in ( ^9|) is independent of cr to write it as ^^(p, g) = ^5(7). By (||), 

(a^ + p2)5^^«-3 ^ („_3)pp"-3^ 

Hence the first term in ( ^9| ) gives, upon integrating by parts over p, 
dp / dq {a'+p^)p^-H'{p - e) 0(7) /«(7) 

dq p"-3 [(n - 3) 6 /« + 2e0 /« + (a^ + p^)[^' /« + /« )] 
dq p"-3 [(n - 1) e /« + (a^ + e^) (0' /« + /« )] , 
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where — dpf*{p, q) and the integrand is to be evaluated at p = e. Inserting the other terms in ( pQI ) and 
simpHfying, we obtain 



{S,,f) = — c.„_i / dq p--' [07' - /• ] • (42) 



Proposition 2 For all n > 3, the regularized source distribution 5^{x + iy) is supported on the oblate 
spheroid x e Eg, and its action on a test function f{x) = f^{'^,cr) is given by 

(43) 

where Ig is the linear functional defined by 



{a^ + e^Y+' r F>{e + iq) 



a^-'^An J^a (e + iq) 



—X T dq, V 



n — 3 



with 



2\v 



F\^)^{a'-q') 



n-2 



(44) 



(45) 



Proof: This follows immediately from ([42|), using 



(a2+e2)(„2_^2) 



Let us verify that as y — > and the source disk shrinks to a point, the source of the singular potential 
contracts to the usual point source. Letting g = in (Esh, we have 



Therefore 



and 



where 



Thus by (g|), 



F«(7) = a"-3(l -e^)'' 



Pil) 



lim4(/«) = 

a— >0 



7/1(7) 
n- 2 

_i (e + M^"-! 



where 7 = e + ia^. 



iPpil) 
n-2 



d^ 



ePpie) 



n-2 



Kn 

A ' 



^ 2 , 



[i-erd^ = B{i^)= 



efHe) 

lim(4,/)=/"(e) + ^A/ 



— An ■ 



Now let e — > 0, and note that since Eg, contracts to the origin, 

lim/''(£) lim / /«(e,cr)dcr = /(O), hence lim lim ( 4, / ) = /(O 



If we assume that the two limits can be exchanged (as will be verified later) and use the definition (|3j) of (5, 
this states that lima^o i^^f) = /(O); giving the following important result. 

Theorem 1 The singular source distribution S{z) is an extension of the usual point-source distribution in 
M", in the sense that 



y ^ S{x + iy) ^ S{x). 



(46) 
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4 Singular Source Distributions 

We are ready at last compute the singular source distributions. By ( ^7| ) and (^), they given by the limit 

(5,/) = lim4(/»). (47) 



However, we cannot simply let e = in the expression (44) for /£(/"), since the resulting integral diverges. 
Recall the decomposition of given in Equation (nj 



Ea^E^UE-UB, (48) 

where iJ* are the interiors of the front and back sides of Eg and B is its rim, the branching sphere. This 
decomposition recognizes the nature of Eg as a limit of closed ellipsoids. Although the two open disks Eg 
look identical to a continuous function, they look distinct to a generalized function like S which is singular 
across Eg. Furthermore, the oblate spheroidal coordinates are an ideal tool for resolving this decomposition 
since g > on g < on E^ and q ^ on B. We will compute {S,f) by decomposing the integral Idf') 
in a way similar to ( ^ ) and then taking the limit e ^ 0. The integral over E^ U E^ gives a sum of single 
and double layer distibutions of dimension n — 1 over the interior of the disk Eg, while the integral over B 
gives a boundary distribution of dimension n — 2 over B. All these distributions are well-defined, giving a 
finite expression for {6, f). 

The divergence of loif') is therefore caused by the attempt to represent the boundary distribution on B as 
part of an integral of dimension n — 1 over Eg, and the above regularization simply amounts to recognizing 
this fact. 



To regularize the integral 



define the Taylor coefficients 



a" ^An J-a (£ + 



r™(e) = —^dl^F'ie + iq) I , = T^{Q) (50) 



and the Taylor polynomials approximating F^{e + iq) and F'{iq) to order g"~^ 



ri-2 ri-2 
m— 771—0 

We rewrite the integral in as 
where 

An.)=/^(^, 0<™<.. (53) 

If /(a;) is differentiable to order n — 2, so is F*[e + iq) = F*{e, q). Then 

F\e + ^q)^Fl_,{e,q)^0{q^^'') 

and the integral on the right-hand side of ( p2| ) has a finite limit as e 0. It therefore remains to compute 
the limit of the sum. We begin by finding 

r7Tz^ = '^r^^ = 2'^"''f ="-2tan-i£^a(e). (54) 
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All the other integrals (^3|) can be computed from the recursion relations 

1 



\m+l 



° {e + iq- e)q"'dq 



{e + iq)'' 

We are interested in the limit £ — > 0, where these relations imply 



Ar(0) = A«_„(0) = A,_„,, 0<m<k. 
The Taylor coefHcients Lk of a{e) are obtained from the expansion 



a{e) 



2tan"^ - 
a 



(-1)' e 



I ^21-1 



Thus 



Ai = TT, A2; — 



2(-l)'+i 
(2^- l)a2i-i' 



^ 2?- Ia2'-i' 



A2;+i =0, ; > 1. 



(55) 
(56) 



(57) 



(58) 



By (|56[), this gives finite values for all the coefficients A™(0). Note that the original expression ( |53| ) diverges 
if we set e = 0. The finiteness of the limits depends on delicate cancellations of contributions from positive 
and negative values of q when e > 0, just as happens in Cauchy's principal value integral. In fact, ( |5^ ) shows 
that the present regularization reduces to the principal value integral when k = 1. 

Using (15^), we can now compute the limit e ^ in (E^): 



= i4(/) + ^ 5]*-™T„,A„_„,_i 

m—O 

r F*{iq)~Fl_^{q) 

An 



Vnif) 



dq. 



(59) 
(60) 



Vn{f) will be shown to be a bounded linear functional of /, and this establishes (5 as a well-defined generalized 
function. 

Since /"(O) is the mean of / over B, the terms in the sum in ( |59| ) represent means of / and its derivatives over 
B. On the other hand, ( ^o|) represents an integral of / and its normal derivative over the interior £'+ U , 
since the boundary terms have already been subtracted in the form of FI^_2{q)- Thus, (|59| ) is the promised 
decomposition of the source distribution corresponding to (|48|). 

Equation (p9|) can be greatly simplified because of certain symmetries satisfied by f'{iq). We claim that 



fi-zq)^P{zq), /« (-zg) = (zg), P^(-^q) ^ -f^{^q). 



(61) 



To see this, note that the coordinates 7 = ztiq denote the same point of Eg, regarded as belonging to Eg. 
Since the test function f{x) is continuous across Eg, so is its integral f^{iq), and this proves the first relation. 
On the other hand, p increases in the ±C direction on E^, which proves the second relation. Finally, the 
third relation follows from the first by differentiation. 

By (|6ll), the function 



F«(zg) = (a2-g2)^ 



(62) 
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is even, so its odd Taylor coefficients vanish: 

m is odd ^ T„, = 0. (63) 



Furthermore, (60) shows that 

n is even ^ Vn{f) = 0. (64) 

Considering the cases of odd and even n separately and inserting the values of from ( |58| ) gives the following 
result. 

Theorem 2 The singular source distribution S(z) is a bounded linear functional which acts on a test function 
f e C""2(]R") as follows. For even n = 2k + 2>A, 

{lf)^^i-l)'T2,, (65) 



where An and T2k are given by ( [45| ) and (50). For odd n = 2fc + 3 > 3, 
where Vn{f) is given by (pQ). 



Proof: To prove (|65[), note first that Vn{f) vanishes by (|64|). Furthermore, the sum in ( p9| ) reduces to a 
single term because TmXn-m-i = unless m is even, which implies that n — m — 1 is odd; but the only 
nonvanishing coefficient A; with odd / is Ai = vr. 



Equation (|66|) follows directly from (|59|), ( |38| ) and ([631). To establish 5 as a distribution, we must still prove 
that Vn{f) is a bounded linear functional for odd n. Since the integrand in ( |60| ) is even when n is odd, we 
have 



T/^n 2»l-"a r F»{zq)-Fl_^{q) ^ 
Vnif) = ^ / — dq. 

An Jo q 



By Taylor's theorem, 



Therefore 



F\^q) - Fl_,{q) = f\q w)-' iF»)(-'\^w) dw. 

(n- 2)! Jq 



n-i.. max|(F«)("-i)(zHI 



\F\iq)^Fl_^{q)\<q-~^M, M = 



(n-2)! 



and 



|14(/)| <2a2M. 

Since M depends boundedly on /, this shows that Ki(/) is a bounded linear functional as claimed. I 

Thus, for even n > 4, the integral over the interior of vanishes and the singular source is concentrated 
on B. This is to be expected, since 7^"" has no branch cut and therefore the only singularity occurs on the 
boundary. 

We can now derive a useful expression for the source distribution 5{x + iy) in the cylindrical coordinates 
(p, (■,"■) adapted to y ^0. As a byproduct, it will be seen that for all n > 3, the test function f{x) need 
only be in C'=(R") with k = [^], instead of C"-2(k«) as assumed in Theorem |. 



11 



Our first task is to rewrite the function F*{iq) in (|6^ ) in terms of cylindrical coordinates. Recall that 

pq 



c = 



p = 



X = pa + Qy, cr e 5"" M 



Then the mean of f on = Ep Hq becomes 



/ h{p<T + Cy)dcT = f{pX), 
which is the mean of / on the sphere Sp,Q — {x : \x ~ ( y \ — p}. 



(67) 
(68) 

(69) 



It will be shown elsewhere [KOO that the source distributions for odd n can be derived from those for even 
n by integrating over one of the coordinates. Specifically, we have: 

Theorem 3 Let n > 3, and denote points in C"+^ by 

z = x + iy = {z, s + it), zeC^. 
For fixed y ^ 0, the singular source distributions Sn{x + iy) in M" and (5„+i(x + iy) in M"+^ are related by 

/oo 
dsSn+iiz,s). (70) 

That is, for a test function fix) in M" we have 

{SnJ)= / ds dx 6n+i{x + iy,s)f{x). (71) 



For y = 0, ( ^0| ) follows from Sn+i{x, s) = Sn{x) S{s). For y ^ 0, this tensor product decomposition fails but 
Theorem || still holds. We therefore state the following theorem for the simpler case of even n > 4, though 
the transformation to cylindrical coordinates derived below will be used later in the explicit computation of 
the extended source in M^. 

Theorem 4 For n = 2k + 2 > 4, the source distribution d is a bounded linear functional on test fucntions 
f{x) in C'^(]R") whose action is given in the cylindrical coordinates {p,(^,cr) by 



T{k 



—,D'pF{p) 



p—a 



where 



" a(p2) 2pdp 

Proof: By (|67|), we have 



^9 7^/ N 2fe-l 



f{p,0) + ^^^-:^fdP,0) 



fp(p^ i) = -zrhi fpip^ C) + - fdp. 0- 



a"^ + p 



(72) 



(73) 



(74) 



Since 



p = ^ q^ = - p^, 
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it follows from that 



F{iq) = {a" - g2 



/«(0,g)+* 



^n — 3 



/(p,0) + z- 



g/^(o,g) 

n - 2 



■/c(p,o) 



(75) 



Note that (ra) makes sense only because F{iq) is even and so does not depend on the sign of g, which cannot 



be recovered from p when C = since sgnq = sgn^. With 

d 1 d 



d{q^) 2q dq 



we have 

oo 

For n = 2A: + 2 > 4, 



m— 



m! ^ 



= —-Dp when p = 0, 
_ (-1)™ 



9=0 ml 



D^F{p) 



27r'=+i 27r 



fc+i 



and (Ba) gives 



-KaUn-l 



1 (^K + 2-* 



Now F{p) already contains one derivative (/^) and computes k more, hence it suffices to have / € 
C'=+i(M"). But the highest derivative of / occuring in ( [72| ) is 

which vanishes at p = a. Hence the highest nonvanishing derivative is of order k, so it suffices for f{p, C) to 
be k times continuously differentiable as claimed, and so does f{x). I 

An application of ( [7l| ) now shows that for n = 2fc + 1, we need / e C'^(M") in order for {6n,,f) to make 
sense. Thus for any n > 3, f needs to be in C'^(R") with 



k 



'n- 1' 




2 





2 ' 



n even 
n odd. 



5 Computation of Sources in R3 and R4 

We now compute the singular source distribution in M"^ explicitly and interpret the result. Equation ( |6^ ) 
becomes 



F\iq) = p{iq) + iqfl,{iq), 
hence Tq = /«(0), Ti = 0, and (|l|) becomes 

F[{q) = p{0). 

Equation ( |66| ) therefore gives 



(76) 



Using the symmetry of the integrand, we obtain the following result. 
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Proposition 3 The singular source distribution 

1 . 
S{x + iy) = -A- — , x + iyeC, 

4:TTJ 

is a bounded linear functional whose action on a test function f{x) = f^{'^,cr) in C^(MP), is given in oblate 
spheroidal and cylindrical coordinates by 

ilk) =Lo + Li+iL2, (77) 

where 

io = /«(0)-/(a,0) 

Li = --I dq = -aj^ ^-^-^^pdp 







L2 = -a dq = - —^===pdp. 

Jo Q Jo Va - P 

Proof: The action in oblate spheroidal coordinates follows immediately from (|66|). To obtain the action in 
cylindrical coordinates, use (|7j). I 

Note that since 

lim f{a, 0) = /(O) and lim Li — lim L2 = 0, 

a — '0 a — >0 a — >0 

( [TTI ) shows that 

y ^ => S{x + iy) ^ S{x). 

This was already seen in Theorem ^ for all n > 3, but that proof was less rigorous because it depended on 
the assumption that the order of the limits a —^ and e ^ can be exchanged. 

We now state some other interesting properties of the expression ( |77| ) which will help interpret its three 
terms. 

Proposition 4 The monopole and dipole moments of S{x + iy) in MP are 

Q= S{x + iy)dx = l, P= / x6{x + iy)dx = —iy. 

Given a point source with general complex coordinates Zs — Xs + iy^ G C'^, the centroid of its charge 
distribution is 

C{zs)= I x~5{x- Zs)dx ^ Zs . (78) 

Proof: To compute Q, apply Proposition |^ with f{x) = 1. To find P, apply it to the vector-valued test 
function 

f{x) = X = pa + Cy- 

Finally, 

C{zs) — Xs S{x ~ Xs ~ iys) dx + {x — Xs) S{x — Xg — iy^) dx = Xs + iy. 
Equation ( [78| ) is a natural extension to C'^ of the formula 

/ X 5{x — Xs) dx = Xs^ Xs<EM?. 
Jk.3 

Proposition || sheds some light on the nature of the source distribution 5. 



z<. 
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• I/O is the mean of / over the rhii B. Since Li = L2 = when f{x) = 1, we see that the "charge" Q 
resides entirely on this rim. 

• Li is an integral of / over Eg which does not involve its derivatives, thus representing a single layer 
distribution on Eg . Actually, since the contributions from the rim are subtracted, the single layer 
resides on the interior of Eg . 

• L2 is an integral of the normal derivative of / over Eq , so it represents a double layer distribution 
which may be regarded as equal and opposite charge distributions on E^ and E^ . This is confirmed 
by P = -iy. 

Proposition 5 The source distribution in C^, 

^(^) = -^W' 

acts on a test function f E C^(R^) in cylindrical coordinates as follows: 

(J,/) =/(a,0)+a/p(a,0)-za/c(a,0). (79) 



Proof: From rt73) we find 



F{p) 



and becomes {S,f) = dpF{a), which gives d? 

As a — > 0, this gives {S,f) /(O), which confirms that S{x + iy) S{x). This property and Q — 1,P 
-iy, C{zs) = Zs , as in Proposition ^, can be shown to hold for all n [ KOC( |. 



The holomorphic potential in C was derived from a different point of view in |K94, Section 11.2], where it 



was shown to decompose into causal (retarded) and anticausal (advanced) "phy sical w avelets." These are 



closely related to the complex-source pulsed beams in the engineering literature [ HF89 |. 



6 Connection to Spacetime and Wave Equations 

It was proposed in Section 1 that point sources in physics be replaced by extended sources based on a 
continuation of physics to complex spacetime. The distribution S{z) in complex space (C") seems like a 
promising model, since it is supported on the disk Eq (when n is odd) or the sphere B (when n is even). 
But to do physics we need to add time to this complex-space formalism. We will first give an argument 
suggesting that time is already included in C"+^, then justfy this by proving that with this interpretation, 
the extended source distribution ^(z) in C"+^ acts as a propagator in spacetime, generating solutions of the 
Cauchy problem for the wave equation from their initial values. 

Let us write 

z = x + iy={z,s + it)eC"+\ zeC"- 
K^{x,s), y = (y,t)eR"+\ \x\ ^ r. 

Analytic potential theory in C"+^ is based on the complex distance function 7(z) defined by 

= z^ ^ z^ + {s + it)^, Re7>0, 

which reduces to the Euclidean metric in R"+^ when y 0, 

2 2 2,2 
7 ^ x = r + s , 
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and to the Lorentzian metric in 



when 



and y ^ 0, 



7 



This suggests that the imaginary part t of the complex space coordinate Zn+i = s + it should be interpreted 
as time. For this reason, we refer to its real part s as the Euclidean time. 



The idea of time as an imaginary space coordinate dates back to Minkowski [ M23 |, who reahzed in 1908 that 
Einstein's new special relativity theory can be based on a unified four-dimensional spacetime. More recently, 
complex spacetime has be come an important concept in quantum field theory S W64 , GJ87 |, twistor theory 
[PR86| and string theory |P98|. Even so, it is generally regarded as a useful mathematical tool rather than 
a fundamental aspect of physical reality. For example, see the discussion in |MTW73, p. 51]. 



In previous work |K77, K78, K87, K90, K94|, I have attempted to make complex spacetime "concrete" by 
giving detailed physical interpretations of the imaginary as well as the real coordinates. This is also the 
present motivation for developing holomorphic potential theory, where the gemetrical and physical signifi- 
cance of y has been emphasized. (The physical significance of y is related to the directivity of the physical 
wavelets associated with (j){x + iy) pF89|, |K94|, |KOC(|.) 



Assume for simplicity that the Euclidean world is completely democratic with respect to the real space 
variables x g M" and the Euclidean time s, so that there is no prefered direction in M"+^. Then if y ^ 0, 
we may choose a coordinate system in which 



(0,t), i=|y|>0. 



(a;, s + it) e 



In that case, the cylindrical coordinates in 



are 



C = X • y = s, p= - 



and a vector x S M."''^^ is represented by 



X = rcr + sy. 



where we do not distinguish between <j E M" and (cr, 0) G M"+^. As before, we denote by da the normalized 
surface measure on S"'~^, so that J dcr = 1 and 



fir,s) 



f{r(T + sy) da 



is the mean of f{x, s) over the sphere \x\ = r. 

For simplicity, we assume to begin with that the dimension n of space is odd, so that n -I- 1 is even and ( 72 ) 
can be applied to 5(z). We will later extend our results to even values of n by applying (|7l|). Using (72) 
with the substitutions 



gives the action of (5(x -|- iy) on a test function /(x): 



a t 



Tik+^) 



uD'rFir) 



where 



1 d 
2r or 



/>,0) + i- 



2kt 



fs{r,0) 



(80) 



(81) 



(82) 



To establish a connection with the wave equation, define the function /(z) on C"+^ by the convolution 

/(z)= / 5(z-x')/(x')dx'. (83) 
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This is an extension of / from R"+^ to C"+^, since 

y -> 5{x + iy - X.') ^ S{x - x') 



Now 



hence 



7(-z) = 7(z) 



(-Z) 



/(x + zy)^/(x). 



5(-z)=<5(z), 



(84) 



/(x + *y) = / 5(x' - X - *y) fix') dx' = 



(85) 



where /'^ is the test function defined by 

r{x')^fix + x') = f{x + x',s + s'). 

Note that in (^sf), it is ^(x' — iy) rather than S{x' + iy) that acts on f^. We will account for this in Equation 
(p6|) below by letting i ^ —i in (p2|). Thus, since z = x + iy = (a;, s + it), (|8l|) gives 



/(a;,s + ii) 

with 

Therefore 



r(fc + i) 



2fe-l 



r(r,0)+* 



2fci 



/?(^0) 



(86) 



r(r^ + s'y)d*T 



/(x + rcr + s'y) dcr = 



/(a; + rcr, s + s') dcr. 



S„-i 



fix 



rcr, s acr 



is the mean of / over the sphere of radius r centered at x and orthogonal to y, and 



ff,ir,0)^ds'nr,s') 
is the mean of fs over the same sphere. 



fsix + r(T,.s)da = ffir,0) 



(87) 



(88) 



Theorem 5 Let n = 2k + 1 > 3 and f e C''+2(R"+i). Then fix, s + it) belongs to C'^iR"''^) as a function 
of ix,t), and it is the unique "classical" solution to the following Cauchy problem for the wave equation: 

dl fix, s + it)^ fix, s + it) (89a) 
liin /(a;, s + it) = /(a;, s), lini 9t/(a;, s + it) = i/s(a;, s). (89b) 

Furthermore, if fix, s) can be continued analytically in s, then fix, s + it) is that continuation. 

Proof: The first of the initial conditions is just the extension property (p^), which has already been estab- 
lished. In view of this, the second initial condition can be written as 



lini ids + idt) fix, s + it) = 0, 



(90) 



which is the CaucTiy-ffiemann equation at t = 0. If (|9C|) were to hold in a complex neighborhood U of sq, 
it would imply that fix, s + it) is holomorphic in s + it G U. Since analytic continuation is unique when 
it exists, this proves our claim that fix, s + it) is the analytic continuation of fix, s) whenever the latter 
exists. 
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Without assuming analyticity, we now prove that f{x, s + it) solves the above Cauchy problem. Let 

u{x,t)^ f{x,s + tt), v{x,r) = f^{r,0), w{x,r)=tff{r,0), (91) 
where the dependence on the parameter s is suppressed. Then ( p6| ) becomes 

^2 _ ^2 



u{x, t) 



T{k + \) ' 



,.2 _ £2 

r'^''^^v{x, r) H ; r^^~^w{x, r) 



2kt 



(92) 



According to (pTI), v{x, r) is the mean of / over the sphere of radius r centered at x, hence v{x, 0) = f{x, s). 
Similarly, ( ^8|) states that w{x,r) is the mean of ifs over the same sphere, so w{x,{)) — ifs{x,s). We 
therefore need to show that u solves the following Cauchy problem: 



utt{x,t) = Au{x,t) 
u{x,0) = v{x,0), ut{x,0) = w{x,0). 



(93) 



By the definition of Dr, the left equation in (pq) is 



d 



k-l 



-dt{jg-^] t''-'vix,t)^cndt 



tdtj 



fc-i 



t^^-^v{x,t), 



(94) 



where 



1 



1 -S - • •(2fc - 1) = 1 • 3 - --(n - 2) 



The right equation in (pw is 



d 

" \ rdr 



2k 



w{x, r) 



(95) 



Letting ^ = r^/2, a straightforward computation shows that for any function G(^), 



d 
rdr 



Thus (p^) becomes 



[{r' ~ t^)G{i)\ = d\ [(2^ - t')G(0] = 2kd\-^G{0 + (2$ - t^) d\G{i) 



2k 



d 
rdr 



k-l 



G+ir^-t^) 



G. 



k-l 



tdt 



k-l 



(96) 



The sum of (|9J) and (|9^) is precisely the solution u{x,t) of the initial- value problem ( |93| ) with n = 2fc+ 1, as 
expressed in terms of spherical means. See John | |J55| , Courant and Hilbert |CH62, pp. 699-703] or FoUand 
| F95| , p. 170]. 

That u g C^(M"'^) follows because v{x,r) G C''"'"^(R") in x and ( |9^ ) contains k derivatives, while w{x,r) G 
(jk+ifj^n-^ in a; and ( p^ ) contains A: — 1 derivatives. Finally, uniqueness of the solution m is a general property 
of the Cauchy problem. I 

As mentioned earlier, the above result can be extended to an even number n of space dimensions by applying 
the recursion relation (70) between the singular source distributions in C" and C"^^. In terms of the solution s 
f{x,s + it) to the Cauchy problem, this amounts to using Hadamard's method of descent | H52 , CH62|. 
Consequently, the same formula ( |85|) gives the solution of the Cauchy problem for the wave equation in R"'^ 
for all values of n> 2. 
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The support properties of S now imply some important attributes of waves in R"'^. For z = x + iy G C"+^ 
with y 7^ 0, recall that 



supp (5(x + iy) = 



isiy) foroddn>3 
\Eo{y) for all other n > 2. 



From ( |85| ) we can therefore immediately draw the following conclusions about waves (solutions of the wave 
equation) u{x,t) in n space dimensions: 



For odd n > 3, u{x,t) depends on the values of u{x + v,0) and ut{x + v,Q) only in an arbitrarily 
thin shell containing the sphere \v\ — t. (We need a shell, rather than the sphere itself, because of the 



derivatives appearing in (|9J) and (|9^).) This is the strong form of Huygens' principle |BC87|, which 
states that u depends on the initial data only on the light cone. 

For all other n > 2, u{x,t) depends on the values of u{x + v,0) and Ut{x + v,0) in the past cone 
\v\ < t. This is the principle of causality, which states that no signal (information, energy) can travel 
with speed greater than c = 1. (If we rescale time hy t ct with arbitrary c > 0, then the maximum 
propagation speed is c.) 

We emphasize that although holomorphy was needed in the dehnition of S, it is not necessary for the above 
relation to the wave equation. This distinguishes our results from all similar results in the literature of which 



I am aware, where holomorphy of the Cauchy data is essential. See Garabedian |G64, pp. 191-202] and Ryan 



[R9C, R90a, El96a 



Theorem ^ and its counterpart for even n can be extended in various ways. 

• Clearly it is not necessary to assume that t > 0, since the support of (5(x + iy) is symmetric with 
respect to y ^ — y. When t < 0, f{x, s + it) is a solution of the "final-value problem" in terms of the 
Cauchy data at t — 0. 

• The Cauchy data /(x) need not belong to C''+^{M:"+^). When / is a distribution belonging, say, to 



some Sobolev space [F95, Chapter 6], then f{x,s + it) is a distributional solution and the derivatives 



in (94) and (B6[) must be interpreted as distributional derivatives. 



To solve the inhomogeneous wave equation 



utt{x, t) - Au{x, t) = j{x, t) 
u{x,Q)—v{x), ut{x,0) = w{x), 



one can apply Duhamel's principle to solutions of the homogeneous equation [F95|. This involves 
integration of S{z') on the truncated solid light cone with < t' < t. However, we will see that 
in the Clifford setting, the current formalism leads directly to solutions of inhomogeneous hyperbolic 
equations, where the time-dependent source is determined by the given function /(x) in Euclidean 
spacetime. 



7 Extension to Clifford Analysis 



Clifford and quaternionic analyses | BDS82 , R9£, GS97, R9S | are generalizations to R" of one-dimensiona l 



complex a nalysis that are proving to be a unifying and very powerful tool in physics [ H66 , KS96 MM98 , 
098 , B99 |. We now show that all the above constructions generalize naturally to this setting. Let C£„ be 
the complex Clifford algebra generated by elements ei, • • • , e„ satisfying the anticommutation relations 



SkSi + ciCk = 2(5, 



kl 



l< kj<n. 



(97) 
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As a complex vector space, C£n has dimension 2" with basis vectors 

Bk = e^^- ■ -e^^ , K = {ki,- ■ ■ ,kp}, < p < n, 1 < ki < ■ ■ ■ < kp < n, 

where the element labeled by the empty set K = (p = 0) is by definition 60 = 1. Thus a general vector 
has the form v — J^k^k^k, where Ck G C and the sum runs over the 2" sets K as above. The element 
J2k=i ^k^k G C'in is identified with the vector z G C", and by ( ^ ) the products and squares of such vectors 
satisfy 

n 

zw + wz = 2 ZkWk =2zw, z"^ = z ■ z = 7(2)^, (98) 
fc=i 

where 7(2;) is the complex Euclidean distance function. This connection will be our basis for generalizing 
holomorphic potential theory to the Clifford setting. 

Clifford analysis is a noncommutative calculus dealing with Clifford-valued functions 
/:M"^a„, f{x)^J2^KfK{x) where : M" ^ C. 

K 

The primary tool is the Dirac operator D = '^^.Skdk, which is closely related to the exterior derivative 



[AM7S] but in addition incorporates the underlying metric. It acts on Clifford- valued functions from either 



left or right by 

k=l k=l 
K 

Thus £) is a "square root" of the Laplacian A in M". It is an elliptic operator because the relations jo^ ) 
are based on the Euclidean metric in K". In 1928, Dirac formulated a similar operator in the Minkowskian 
spacetime M'^'^, where it is hyperbolic and its square is the wave operator. This formed the basis for his 



relativistic wave equation of the electron |D58|, which had a revolutionary impact on physics, including 
especially the dramatic prediction of antimatter. Mathematicians usually prefer Euclidean Dirac operators 
because, among other things, they yield powerful methods for solving boundary-value problems generalizing 



those in one-dimensional complex analysis |GS97|. For this and similar reasons, most mathematical work 



on Dirac operators is restricted to the elliptic case; see the discussion in |098, p. 1]. Consequently, rigorous 
analyses of Maxwell's equations by Clifford methods usually deal with static or time-harmonic fields |KS96 
MM9^ . By assuming that the boundary/Cauchy data is holomorphic, it is possible to arrive at solutions 



of hyperbolic Dirac equations through analytic continuation [ |R90| , R90a, El96a|, generalizing the method 
employed by Garabedian for establishing a connection between solutions of the Laplace and wave equations 
|G64|. As in Section the present method is not restricted by this assumption. We will "Cliffordize" 



the holomorphic potential theories in C" and C"^^, interpreted respectively as complex space and complex 
spacetime. The relation established in Section ^ between Laplacians and wave operators then yields the 
desired connection between elliptic and hyperbolic Dirac operators. 

The Clifford counterpart of the Newtonian potential (f>{x) is the Cauchy kernel C : M" C£n , defined as 
the fundamental solution of D: 

DC{x) = C{x)D ^6{x), limC(a;) = 0. 

r — >OQ 

Because = A, the solution is easily expressed in terms of the Newtonian potential: 

C'ix) = Dcj^ix) - (101) 



w„(2 - n) 
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The expression on the right remains vaHd ioi n — 2 if we take 4> — (27r) ^ Inr. Replacing r{x) with 7(2;) 
immediately gives an extension to C", where the point source S becomes the extended source 6: 

C{z) = D4>{z) = (t){z)D = , 2;eC", n>2 (102) 

^ DC{z) = C{z)D = A(l){z) ^ S{z). (103) 

For all even n > 2, C{z) is holomorphic on the complement of the null cone J\f. For odd n > 3, it inherits 
the branch cut from 7. Thus for all n >2, C{z) is holomorphic on the complement of the set 



{x + iy \ {x + iyY = 0} for even n > 2 
{x + iy e C" I X e Ea{y)} for odd n > 3. 



Given a Clifford-valued test function / : M" ^ C£„ , we define its extension to / : C" — > C¥„ exactly as 
before: 

f{z)= [ ~d{x' -z)f{x')dx', ^ lim fix + iy)^ fix). (104) 



Recall that for / to be defined, / must be C'' with k = [^V'"] n> 3. It can also be shown | KOC | that / 
must be if n = 2. Inserting one of the expressions ( |103| ) for 6 and integrating by parts gives 

fiz)=J Cix' - z)Dfix')dx' ^ - j Cix' - z)Dfix')dx\ (105) 

where D denotes the left-acting Dirac opeator with respect to x' . We will use this expression to derive an 
extended version of the Borel-Pompeiu formula. Let M be a bounded domain in R" with piecewise smooth 
(C^) boundary dM, and let 



fM{x) = Xm{x) fix), where XMix) 



1 ifxe M 

if a; e M' = K"\M. 



We do not need to define on dM. Taking the distributional derivatives of /^j gives 

Df,,ix) = iDxMix))f{x) + XMix)Dfix). (106) 



We want to substitute this into (105) to obtain an expression for the extension f^iiz) of fM{x). This will 
make sense if the singularities of Df^ do not meet the singularities of C{x — z). That will be the case if we 
assume that 

z i dM + Sn = {xb + z\xbe dM, z e 5„}, (107) 
and we refer to such points z G C" as regular with respect to dM. Note that when z E R", this means 



simply that z ^ dM. For z regular, we may substitute ( |lOq ) into (105) to obtain 

fMiz) = - f Cix' -z)iDxM(x'))f{x')dx' - f Cix' ^z)xM{x')Dfix')dx'. (108) 

JR" JR" 

We claim that the first term can be written in classical (non-distributional) form as 

Cix' ~ z)nix')fix')daix'), (109) 

dM 

where nix) is the outgoing unit normal at a; G dM and daix) is the area measure on dM induced from the 
volume measure dx. To see this, note that there exists a differentiable function /i(a;) such that 

M = {x I ^ix) > 0}, dM = {x I nix) ^ 0}, V/i(a;) = -\\/nix)\nix). 
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Then Xm can be expressed in terms of the Heaviside step function 9 by Xm{x) = 0{fi{x)), hence 
DxMix) - Si^iix))D^i{x) = S{^i{x))V^i{x) = -S{^i{x)) \V fi{x)\ n{x) . 



The connection to the expression (109) can now be made by using the imphcit function theorem with fJ.{x) 
as one of the local coordinates. The second term in ( |108| ) reduces to the integral over Af , giving the following 
result. 



Theorem 6 (Extended Borel-Pompeiu formula) 

Let z e C" be regular with respect to dM . Then 



[ C{x' - z)n{x')f{x')daix')- [ C{x' - z) Df{x') dx' . 

JdM JM 



(110) 



For z x IE M"\9Af , Equation ( |l04| ) applied to /m reproduces fM{x) since S{ x' — x ) = 5{x' — x). Thus 
(110) reduces to the usual Borel-Ponipeiu formula for Clifford- valued functions |GS97|, 



[ C{x' -x)n{x')f{x')da{x')- [ C{x' - x) Df{x') dx' = \ 

JdM J M y 

We may interpret f{x) as a field generated by the source function 

3{x) = Dfix). 



f{x), xeM 

0, xe M'. 



(Ill) 



(112) 



Then (111) solves the boundary- value problem for ( |ll2| ), expressin g th e field / inside M in terms of its 
sources in M and its value s on d M. (To investigate the limit of (111) as a; ^ dM, one also needs the 
Plemelj-Sokhotzki formula GS97 |.) In particular, / is said to be left-monogenic in M ii j{x) = in M. In 
that case, the second term in (111) vanishes and the Borel-Pompeiu formula reduces to a multidimensional 
generalization of Cauchy's integral formula, with monogenicity replacing holomorphy. 

To describe waves, such as a time-dependent electromagnetic field, we ascend to C""*"^ as explained earlier. 
In the notation of Section ^, let 

z = x + iy eC", z = (2,zo) = x + iyeC"+^ 
y — {0,t) x—{x,s) z = {x,s + it). 

The generators of Cln+i are {eo, • • • , e„}, satisfying relations identical to ( p7|) but with < fc, Z < n. The 
elliptic Dirac operators in R" and R"+^ are 

n n 

D = ^ekdk, B = ^ekdk ^ D + eods , do = ds ■ 

k=l k=Q 

We also define the hyperbolic (space-time) Dirac operator by formally substituting s ^ it in D: 

D = £> — ieodt ■ 

Then 



A -I- a? and 



where A is the spatial Laplacian in M" and □ is the wave operator in M"'^. 

Theorem 7 (Inhomogeneous Dirac/Msixwell Equation) For n = 2fc + 1 > 3, let f : 

be a (7'^+^ function and define the functions j : M"^"'^ — > C£n+i and j : C"^^ — > C£n+i by 

^f{x, s) = j{x, s), t)f{x, s + it) ^ j{x, s + it). 



C£n+1 

(113) 
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Then j{x, s + it) is a solution of 

t)j{x, s + it) ^ 0, 
and the system {f,j} satisfies the initial conditions 

f{x,s) = f{x,s), j{x,s) = j{x,s). 



(114) 



(115) 



Proof: This follows directly from Theorem ^. Since f{x,s + it) is a solution of the homogeneous wave 
equation, j{x, s + it) is and 

Dj = &f = □/ - 0. 
From the initial conditions on /, we have f{x, s) ~ f{x, s) and 

\imj{x, s + it) = Df{x, s) — ieo lirn dtf{x, s + it) 

Df{x, s) + eodsfix, s) = D/(a;, .s) = j{x, s). I 

For n = 3, the right side of (|113| ) is precisely the Clifford form of the inhomogeneous, time-dependent 
Maxwell equations [H66, 098, 699] ], where 



f{x,s + it)=^ ^ efj,^F^^{x,t) = F{x,t) 



(116) 



is the electromagnetic field (a bivector). 



j{x, s + it) ef,Jf,{x, t) = J {x, t) 

fj.=0 



(117) 



is the charge-current density (a four- vector) , and the Euclidean time s is a free parameter. As with our 
distance function, the Minkowskian metric appears in (116) and (117) when the time components Jq and 
Fok {k = 1,2,3) are imaginary and the spatial components are real. 

Again, we did not need to assume that / or j are holomorphic in s + it. Note that although j{x, s + it) is a 
time-evolved extension of j{x, s), it can be shown |KO0| that it is not the one obtained by convolving with 
S. Instead, the time dependence of the charge-current density is governed by (114), the scalar part of which 
is the continuity equation implying conservation of the total charge. 



In (|TT|) and (|lj), we have taken M = M"+i for simplicity, giving a pure initial-value problem. These 
equations may also be formulated in a bounded region M C R""''^ by using the extended Borel-Pompeiu 
formula for /^(a;, s + it). This suggests that Theorem ^may be used to solve mixed initial/boundary value 
problems, provided one is careful about dealing with the singular points {x,s + it) g dM + Sn+i ■ These 
are in the domain of influence of the boundary; that is, they can be reached at time t by signals originating 
from the boundary dM at time t — 0. However, this still does not give a mechanism for boundary effects 
initiated at t > (such as reflections) to influence the solution. This and related questions will be treated 
elsewhere. 
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